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Abstract 



We point out that the decay mode D° — > K°K° has no factorizable contribution. In the chiral 
perturbation language, treating D° as heavy, the 0(p) contribution is zero. We calculate the 
nonfactorizable chiral loop contributions of order 0(p 3 ). Then, we use a heavy-light type chiral 
quark model to calculate nonfactorizable tree level terms, also of order 0(p 3 ), proportional to 
the gluon condensate. We find that both the chiral loops and the gluon condensate contributions 
CNJ ' are of the same order of magnitude as the experimental amplitude. 
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^ ■ 1 Introduction 

a: 

O ,| The decay mechanism of the weak nonleptonic D° decays has motivated numerous studies [|]]- 

[pTofl - For nonleptonic decays of D mesons, as well as for ICs and -B's, the so called factorization 
hypothesis has been commonly used. For nonleptonic decays, the effective Lagrangian at quark 
level has the form 



-IF 



^aQ,, (i) 



where the coefficients Cj contain all the short distance electroweak and QCD effects to a certain 
order in perturbation theory, and the Qj's are quark operators. Typically, these quark operators 
are products of (pseudo) scalar- or vector- currents: Q = j{2). Then, for a nonleptonic decay 
M — > Mi + M2 , the factorization hypothesis (-also called vacuum saturation approximation) gives 
prescriptions of the form 

(Mi M 2 \Q\M) - (Mi|j(l)|0)(M 2 |i(2)|M) . (2) 



The factorization hypothesis are known to fail badly for nonleptonic K decays [11, 12, |13|. On 
the other hand, there are certain heavy hadron weak decays where factorization might apply. 
Recently, the understanding of factorization for exclusive nonleptonic decays of B mesons in terms 
of QCD in the heavy quark limit has been considerably improved Q]. In this paper we will discuss 
nonfactorizable terms for D decays, in particular for the decay mode D° — ► K°K°. 
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Even though the factorization hypothesis might work reasonably well if one is interested in an 
order of magnitude estimate, it does not reproduce experimental data completely. For example, a 
naive application of factorization in charm decays leads to rates for the D° —* vr ^ , D° — ► 7r°7r°, 
D° — ► K + K~ D° — > tt + tt~ decays which are too strongly suppressed. Moreover, and this is 
the important point of this paper: in D° — ► K°K°, factorization misses completely, predicting a 
vanishing branching ratio, in contrast with the experimental situation. 

To see this, note that at tree level the D° —* K°K° decay might occur due to two annihilation 
diagrams Q which could potentially create the K°K° state. However, they cancel each other by 
the GIM mechanism. Moreover, in factorization limit, the amplitude is proportional to 

(K^^V^iOlA^D ) ~ (pko-p^^dp 1 ^ = . (3) 

In many of the studies (e.g. g, [3|, [|, ||, [j]]) this decay has been understood as a result of final 
state interactions (FSI). In the analysis of ref. M the rescattering mechanism included K + K~ and 
7r + 7r~ states leading to a branching ratio B(D° — > ET°X°) = ^B(D° — > K + K~). Experimental 



data on the other hand are ||] 5(Z) = (6.5 ± 1.8) x 1(T 4 and 5(£>° -> K + K~) = 

(4.25 ± 0.16) x 1CP 3 . A recent investigation of the D° — » K°K° decay mode performed in || has 
focused on the s channel and the t channel one particle exchange contributions. The s channel 
contribution has been taken into account through the poorly known scalar meson /q(1710) and was 
found to be very small, while the one particle i-exchanges yielded higher contributions, with pion 
exchange being the highest. In the approach of || the D° -> K°K° was realized through the scalar 
glueball or glue-rich scalar meson. 

The second instructive case concerning the factorization hypothesis, is offered by the analyses 
of nonleptonic K meson decays. Namely, it is well known that factorization does not work in 
nonleptonic K decays. Among many approaches the Chiral Quark Model (xQM) |l(| was shown 
to be able to accommodate the intriguing AI = 1/2 rule in K — * tttt decays, as well as CP violating 
parameters, by systematic involvement of the soft gluon emission forming gluon condensates and 
chiral loops at 0(p 4 ) order |l2|]. In the xQM [17|, the light quarks (u,d,s) couple to the would- 



be Goldstone octet mesons (K, n,rj) in a chiral invariant way, such that all effects are in principle 
calculable in terms of physical quantities and a few model dependent parameters, namely the quark 



condensate, the gluon condensate and the constituent quark mass [12, [L6|, pj|]. Also in "generalized 



factorization", it was shown [13| that the inclusion of gluon condensates is important in order to 
understand the AI = 1/2 rule for K — ► 2ir decays. 

As the xQM approach successfully indicated the main mechanisms in K — > tttt decays, it seems 
worthwhile to investigate decays of charm mesons within a similar framework. However, in the case 
of D meson decays one has to extend the ideas of the xQM to the sector involving a heavy quark 
(c) using the chiral symmetry of light degrees of freedom as well as heavy quark symmetry and 
Heavy Quark Effective Field Theory (HQEFT). Such ideas have already been presented in previous 
papers [19, |2^, ^] and lead to the formulation of Heavy-Light Chiral Quark Models (HLxQM). In 
our formulation of the HLxQM Lagrangian, an unknown coupling constant appears in the term 
that couples the heavy meson to a heavy and a light quark. Our strategy is to relate expressions 
involving this coupling to physical quantities, as it is done within the xQM J12|J . We perform the 
bosonization by integrating out the light and heavy quarks and obtain a heavy quark symmetric 
chiral Lagrangian involving light and heavy mesons [22, 23 1. 



Because the 0{p) (factorizable) contribution is zero as seen in Eq. (|3|), we will try in this paper 
to approach to the D° — > K°K° decay systematically to 0(p 3 ). We do this by including first the 
nonfactorizable contributions coming from the chiral loops. These are based on the weak Lagrangian 
corresponding to the factorizable 0(p) terms for D — > tt + tt~ - and D° -> K+K~. (Note that the 
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velocity = Pp/mc is considered to be 0(p ) in the chiral counting). Second, we consider the 
gluon condensate contributions, also of 0(p 3 ) within the xQM and HL^QM framework. It should 
be noted that because the energy release in D — > KK is of order one GeV (-in contrast to 200 
MeV for K — > 2ir), the next to leading 0(p 5 ) terms might be almost of the same size numerically 
compared to our 0(p 3 ) terms. Still, the amplitude of D° — > K°K° calculated within the framework 
of 0(p 3 ), has a reliable order of magnitude. Note that we have also omitted I/tuq terms in the 
framework of HQEFT. 

Our paper is organized as follows: In Section 2 we write down the basic Lagrangians including 
the weak Lagrangian at quark level (with special emphasis on the terms giving rise to the nonfac- 
torizable gluon condensate contributions) as well as the standard strong chiral Lagrangians for the 
light and heavy meson sectors. The chiral loop contributions to the decay amplitudes are presented 
in Section 3. The details of the Heavy - Light Chiral Quark Model (HLxQM) are presented in 
Section 4, while the bosonization of the weak quark currents is given in Section 5. The results are 
given in section 6. Appendix A contains some details from the chiral loop integrals, Appendix B 
some details about the D meson decay constant, while Appendix C contains some loop integrals 
within the HLxQM. 

2 Basic Lagrangians 

The effective weak Lagrangian at quark level relevant for D — > -kit, KK is 

C W = G [c A (Qa - Qc) + c B (Q^ ~ Q ( b)} > ( 4 ) 
where G = — 2y/2G F V us V* s , and 

Qa = (sl^cl) (ul7 m s l ) , Q c = (d L ^c L ) (uL7^d L ) , 

Q ( b ] = («lVcl) (q^QL) , (q = s,d), (5) 

are quark operators. 

Using Fierz transformations and the following relation between the generators of SU(3) C (i, j, I, n 
are color indices running from 1 to 3 and a is an index running over the eight gluon charges): 

SijS ln = l-JnSij + 2 C tf 3 , (6) 



one obtains 



Qa = ±-Q { b ] + 4 S) - Qc = ^-Qf + < ] , 



(7) 



Q'b' = Wc Qa + Ra, Qb = jj-Qc + Rc, 
where the i?'s correspond to color exchange between two currents and is genuinely nonfactorizable: 

R A = 2{s Ll n a c L ) (u L t a ^s L ) , R c = 2(d Ll n a c L ) (u L t a ^d L ), 

R { § ] = 2(u Ll n a c L ) (q L t a ^q L ) ; (q = s,d), (8) 

The operators can be written in terms of currents, for instance: 

n (s) _ n (d) _ jY -n „(s) _ R (d) _ jY ,a ./i,o /qn 
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where 



= u L 7 M t a c L ■ jT = 7 M s L - [s ^ d] . (10) 

The currents without color index are given by the corresponding expressions dropping the color 
matrix. 

The factorization approach amounts to writing the currents , jj£ in terms of hadron (in 

our case meson) fields only, so that the operator Q^' — in the left equation of (§) is equal 
to the product of two meson currents. The color currents in @ are then zero if hadronized 
(mesons are color singlet objects). There is also a replacement of the Wilson coefficients in the 
hadronized effective weak Lagrangian ca,b ca,b(1 + 1/N C ). Combining heavy quark symmetry 
and chiral symmetry of the light sector, we can obtain the weak chiral Lagrangian for nonleptonic 
D meson decays due to factorizable terms. Then we can first use this to calculate nonfactorizable 
contributions due to chiral loops. Second, we can calculate the color currents' contribution using 
the gluon condensate within the framework of the HL^QM. 

Treating the light pseudoscalar mesons as pseudo-Goldstone bosons one obtains the usual 0(p 2 ) 
chiral Lagrangian 

42 = y tr(d^d^) + ^ tr(M q E + M 9 tf) , (11) 

where S = exp (2i$//) with $ = ^7 X^ir 3 containing the Goldstone bosons it, K, rj, while the trace 
tr runs over flavor indices and J^Aq — diag(77i u , md,m s ) is the current quark mass matrix. From 
this Lagrangian, we can deduce the light weak current to 0{p) 

J? = -*ytr(£d M £tA x ) , (12) 

corresponding to the quark current = qLl^ X Ql- (A is a SU(3) flavor matrix.) 

In the heavy meson sector interacting with light mesons we have the following lowest order 0(p) 
chiral Lagrangian 

C^. = - Tv(H va iv ■ D ab H vh ) - g Tr ■(3 va H vb 'y li A% a 7s) , (13) 

where iD^H b = id^H a - H b V£ a , the trace Tr runs over Dirac indices. Note that in ( |i~3|) and the 
rest of this section a and b are flavor indices. 

The vector and axial vector fields and A^ in (|l3|) are given by: 

V„ = \{td»t ] + t ] d»t) ; A = -(£8^ - , (14) 

where £ = exp (i$//). The heavy meson field H va contains a spin zero and spin one boson: 

H va = P + (P^-iP 5al5 ) , (15) 



H, 



7 U (^a) T 7 



pt ^ 



5a75 



P 



+ ! 



(16) 



with P± = (1 ±7^f At )/2 being the the projection operators. The field P^(P^) annihilates (creates) 
a pseudoscalar meson with a heavy quark having velocity v, and similar for spin one mesons. 
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For a decaying heavy quark, the weak current is given by 

J* = q al X LQ, (17) 

where L = (1 — 7s)/2 and Q is the heavy quark field in the full theory, in our case a c-quark, and 
q is the light quark field. (For flavor a = it, this is the current in (§).) 

From symmetry grounds, the heavy-light weak current is to O(p ) bosonized in the following 
way [23] 

J^anTr^LHv^l], (18) 
where ajj is related to the physical decay constant fu through the well known matrix element 

<0|n 7 A 75c|£°> = -2(0] J*\D°) = im D v x f D . (19) 
Note that the current fllq) is 0(p°) in the chiral counting. 

3 Chiral loop contributions 

In the factorization limit there are no contributions to D° — > K°K° at tree level. The observation 
of a partial decay width B(D — > /f ^ ) = (6.5 db 1.8) x 10 -4 on the other hand implies that we 
can expect sizable contributions at the one loop level. Calculations to one loop in the framework 
of combined chiral perturbation theory and HQEFT involves a construction of the most general 
effective Lagrangian that has the correct symmetry properties in order to make the renormalization 
work. We discuss constructions of counterterms in the end of Sect. 5. 

We work in the strict MS renormalization scheme, where we put A = - — 7 + ln(47r) + 1 equal 
to one in the loop calculations. This choice, A = 1, determines the appropriate renormalization of 
couplings in the C(p 3 ) effective Lagrangian and is the same as made by Stewart in 28], while it 



differs from the one used by authors of Ref. [22], who use A = 0. We consider only contributions 



coming from the ca part of the weak Lagrangian as cb is suppressed compared to ca [pq ]- 

Writing down the most general one loop graphs with two outgoing Goldstone bosons (K° and 
K°) one arrives at 26 Feynman diagrams. A number of these give zero contributions and are shown 
on Figures [l],§,|3], while the graphs that do contribute to D — ► K°K° decay are shown on Fig. |3|. 
Note that factorizable loops which renormalize vertices are omitted (they do appear, however, in 
the loop determination of the an coupling related to fjj. See Appendix B.) 

To shorten the notation, the common factors in the S matrix have been organized such that 
the amplitude is written 

M(D° -> K°K°) = - ^ c A V us V* > (20) 

where F = ^ n F n is the sum of the amplitudes corresponding to the graphs on Figs, [j], ||. [|. 
The partial decay width for the decay D — ► K°K° is then 

1 r 2 ip|2 

1 do^koko ~ ^^ C A\Vus V cs \ j^p , (21) 
where p is the K°(K°) three-momentum in the D° rest frame 



P=2\l m D-^ m K ■ (22) 
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Figure 1: Diagrams that give zero contribution since the relevant vertices appearing in the heavy 
meson chiral Lagrangian ( |l3|) are zero. The double line represents heavy meson D or D*, while 
dashed lines denote pseudo-Goldstone bosons. 
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Figure 2: Diagrams that give zero contributions since the loop integrals are zero. The double line 
represents heavy meson D or D*, while dashed lines denote pseudo-Goldstone bosons. 
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^ / \ t t 
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Figure 3: Power suppressed diagrams (neglected in the calculation). 



Fx) F 2 ) F 3 ) K+(tt+) 



D° / \ D° / \ D° / \ 



' \ / 

/ S / 

/ I I 

I I I 

D° \ j D° \ 



K°(K°) ' N K°(K°) 
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Figure 4: The nonzero diagrams in D° — > decay. 
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Figure 5: The momenta flowing in the graphs corresponding to a) power suppressed Fg amplitude 
and b) the leading contribution F4 amplitude. 



The nonzero amplitudes corresponding to the graphs on Fig. || are 

ga H 13 
/ 



F 1 + F 2 + F 3 = -^-[A^K.A^) - A*Ji(m K ,A*)], (23) 



F 



■{(m 2 D - 2rn%)[No(m 7r ,mj D ) - N (mK,m 2 D )} + 



(26) 



3/2 

+ ml[N 2 (m«,ml)-N 2 (m K ,m 2 D )}+ (24) 
+ [N 3 (m 7V ,m 2 D ) - N 3 (m K ,m 2 D )] - (ml - m 2 K )Ng(m^,m 2 D )^ , 

F 5 + F 6 = ^^l[/ x ( m7r ) - h {m K )l (25) 
F 7 + F 8 = ^[A d (J 1 (m K ,A d ) + J 2 (m K ,A d )) - A s (j 1 (m 7r , A s ) + J 2 (m w , A s )) 

+ m D ^-I 2 (m K ,A d ) - mD^-him-x, A s ) + -^-/i(m^) - — =-ii(m,r) 
Ad A s 2A^ 2A S > 

where A^ = m (*) —m D o and A„ = rri£)/2 + Ag for q = d,s. Note that A„ are of the order m^/2, a 

consequence of relatively high momenta flowing in the loops of graphs Fj,F 8 . The functions I\(m), 
I 2 (m,A), Ji(m, A), J 2 (m, A), Ng(m,k 2 ), N 2 (m,k 2 ), N 3 (m,k 2 ) appearing in the amplitudes (24-27) 
can be found in Appendix |A]. 

It should be noted that in eq. (23-26) all the expressions vanish in the exact SU (3) limit, where 
rriK — ► and A s — > A d , A s — ► A^. This shows explicitly that the -D — ► K°K° decay mode is a 
manifestation of SU(3) breaking effects (as already noted by H. Lipkin if U symmetry is exact, 
then T(D° -> K°K°) = 0). 

The amplitudes shown on Figs. ||J2|,^ are either exactly zero or are suppressed by powers of 1/mn 
and g = 0.27. The amplitudes corresponding to diagrams on Figs. ||,^ are zero due to symmetry 
reasons (because there are no such couplings in the heavy sector chiral Lagrangian (|l3|), or because 
of Lorentz covariance), while the amplitudes Fg, Fiq and Fu shown on Fig. [|are power suppressed. 
An analysis of the loop integrals leads to the conclusion that Fg ~ g(q/mD) F4, Fiq ~ g{q/friD)F^ 
and -Fn ~ g^ (q/mj^F^, where q is a typical loop momentum less that mfl/2,- so the suppression 
need not be substantial. However, a direct evaluation of the amplitude F10 shows that it is about 
20 times smaller than F4. Therefore, in our numerical calculation we neglect contributions of Fg, 
F\g and Fu. 



4 A Heavy-Light Chiral Quark Model (HLxQM) 

The nonfactorizable contributions to D° — > K°K° coming from the chiral loop correction at the 
meson level obtained in the previous section are not the only contributions to 0(p 3 ). In the effective 
weak Lagrangian @ there are, after Fierz transformations, terms that involve color currents (see 
(p|) , (|i0|) ) . As mesons are color singlet objects, the product of color currents does not contribute 
at meson level in the factorization limit. However, at quark level they do contribute through the 
gluon condensate as will be shown in the next section. In order to estimate this contribution we 
have to establish the connection between the underlying quark-gluon dynamics and the meson level 
picture. This is done through the use of the Heavy-Light Chiral Quark Model (HL^QM). 
Our starting point is the following Lagrangian containing both quark and meson fields: 

£ = C-HQ + £-xQM + £lnt , (27) 

where 

C HQ = Q v iv DQ v + 0{m Q l ) (28) 

is the Lagrangian for Heavy Quark Effective Field Theory (HQEFT). The heavy quark field Q v 
annihilates a heavy quark with velocity v and mass ttlq. is the covariant derivative containing 
the gluon field. The light quark sector is described by the Chiral Quark Model (xQM): 

C xQM = q.{i'fD li - M q )q - m x (q R ^q L + q L Hq R ) , (29) 

where q = (u, d, s) are the light quark fields. The left and right-handed projections q^ and q R are 
transforming under SU(3)l and SU(3)r respectively. M q is the current quark mass matrix, and £ 
is a 3 by 3 matrix containing the (would be) Goldstone octet (jr,K,r)), appearing already in (|TT|). 
The quantity m x is interpreted as the (SU(3)-invariant) constituent quark mass for light quarks, 
supposed to appear due to chiral symmetry breakdown at a scale A x ~ 1 GeV. 
The xQM has a "rotated version" with flavor rotated quark fields \ given by: 

xl = €Ul ; xr = £qr ; £ • £ = s • ( 3 °) 

In the rotated version, the chiral interactions are rotated into the kinetic term while the interac- 
tion term (proportional to m x in (f29|) and responsible for the tt -quark couplings) become a pure 
(constituent) mass term: 

£ X QM = X WiiD^ + V M + 75^) - m x ] X - xM qX , (31) 

and M q defines the rotated version of the current mass term: 

M q = $M q $ R + £Mj£ L = MKR + MfL = + M g A 75 , (32) 

where L = (1 — 7s)/2 is the left-handed projector in Dirac space, and R is the corresponding 



right-handed projector. The Lagrangian (31) is manifestly invariant under the unbroken symmetry 
SU(3)v (if Mq is formally chosen to transform as £). In the light sector, the various pieces of the 
strong chiral Lagrangian (11) can be obtained by integrating out the constituent quark fields x- 



This is the bosonization to be discussed in more detail in the next section. 

Similarly, a left handed current can be written (A^ is a SU(3) flavor matrix) 

q L r\ X q L = XL1^ X XL ; A X = £ f A X £ ■ (33) 
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By coupling the fields A^, M q ,A , A x to quark loops, the chiral Lagrangians of the weak sector can 
be obtained. 

In the heavy-light case, the generalization of the meson-quark interactions in the pure light 
sector xQM is given by the following SU(3)y invariant Lagrangian [19, 2C, 21, [24|]: 



C 



int 



-Gh [x.f H v f Q v + Q v H v f xf] 



(34) 



where Gh is a coupling constant which is related through bosonizations to physical quantities like 
an and g appearing in (|l^) and (|l8|), as well as /„- and m x . (See Appendix C). 

Within HQEFT the heavy-light weak current in ( |T7| ) will, below the renormalization scale 
= m c , be modified in the following way [pi]] : 

Ja = C^) X bti a l X LQ v + C v (fi) X bZi a v X LQ v , (35) 

where the coefficients C~* v are determined by QCD renormalization for [i < m c . However, for 
~ A x , C 7 ~ 1 and C v ~ 0. The bosonization of p5|) will lead to (p^|) by using (|34"|). 



5 Bosonization 



The Lagrangian (^) from the previous section can now be used for bosonization, i.e. we integrate 
out the quark fields. This can be done in the path integral formalism, or as we do here, by expanding 
in terms of Feynman diagrams. For instance, the lowest order (kinetic) chiral Lagrangian (|ll|) in 
the light sector (involving tt, K, ry's) can be obtained by coupling two axial fields to a quark loop 
using the Lagrangian in Eq. (|31 



d d p 
(27r) d 



Tr [(7,75-4") S{p) (l P ^A p ) S[p)} ~ Tr [A^] 



(36) 



where Sip) = (7 -p ■ 



, and the trace is both in flavor and Dirac spaces. This is the standard 



form of the lowest order chiral Lagrangian (|l l[) , which can easily be seen by using the relations 

Aa = -U(d^)£ = l£t(a E) £t. (37) 



Similarly one obtains the lowest order Oip) strong chiral Lagrangian (13) in the heavy sector . 

Let us now consider the bosonization of the pure light weak current. The lowest order term 
Oip) is obtained when the vertex from ( |33"| ) and axial vertex (~ A^) from (31) are combined 
with quark loops (see Fig. |6|): 



iN, 



I Wf Tr i(^ LAX ) S (P) (7*76^) Sip)] ~ Tr [A X A, 



(38) 



This coincides with (II4 ) when (^) is used. 

To obtain a nonzero nonfactorizable contribution to D° — ► K°K° at tree level, we have to 
consider the color current j^' a to 0{p 3 ), involving insertions of the "mass fields" A4 q in (32). From 
Fig. [?], one obtains the contribution: 

d d p 



j*' a iG b ,A,M q \Fig. 7) 



(2vr) d 



Tr 



l^L A ) Sip) iia^A*) Si P )M q S 1 ip,G b ) 



(39) 
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Figure 6: Feynman diagram for bosonization of left-handed current to order 0(p) 




Figure 7: Diagram for bosonozation of the color current to 0(p 3 ) 



where 



S 1 (p,G b 



9s n b Jo 



ad 



(40) 



is the light quark propagator in a gluonic background (to first order in the gluon field) and g s is 
the strong coupling constant. Moreover, a, b are color octet indices. Summing all six diagrams with 
permutated vertices compared to the one in Fig. [?] we obtain in total: 



jf a {G\AM q ) 



I2m x 16?r 2 



G a ' KX [ie^x T x > p + ( WpA 



where (We have used the analytical computer program FORM PQ] ) 
The S's are chiral Lagrangian terms: 



Si 



Wp^Tf'f] , (41) 



(42) 



S L p = Tr 



k X ApM q 



2i 



-Tr 



\ x {D p Y>)M\ 



Sf = Tr 



ft 



2i 



-Tr 



X x M q (D p ^ 



A A [A p M q R + M q L A p 



\ x ((D p E) J$M q Y$ - Y,M\Y,{DpYt 



(43) 



Within the heavy-light sector, the weak current can be bosonized to lowest order (O(p )) by 
calculating the Feynman diagram shown in Fig. ||, left. The obtained result is Eq. (^) with an 
related to Gh (see Appendix C). 
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7^ Lt a 



Figure 8: Diagrams representing bosonization of heavy-light weak current. The boldface line 
represents the heavy quark, the solid line the light quark. 



The bosonization of the color current given by ( p5|) with an extra color matrix t a inserted and 
with an extra gluon emitted is given by the following loop integral (Fig. ||, right): 

d d k 



J a {H v G a ) f 



(27T> 



7 Tr 



{-iG H H v ^) f {iS 1 ( i k,G b ))(j' 7 Lt a ){iA v (k)) 



(44) 



where A v (k) = P+/k ■ v is the heavy quark propagator. Notice that emission of a gluon from the 
heavy quark is suppressed by 1/wiq and omitted. The result can be written 



r{H v G a )f = G H g s G a a0 Tr la L (H v tf)f (l G1 a^ -ilea (7V 



(45) 



where 1q\ and Iq% are loop integrals given in Appendix C. Keeping only the pseudoscalar field P5 
representing D°, we find 



(46) 



where B> e ,g are obtained from loop integrals in ([45]) . Then we find the nonfactorizable (gluon 
condensate) contribution: 



£ e// (L>°decay)( G 2) 



1 



1 16tt 2 J \12m x 16vr 2 



X V, 



X,p 



+ B 9 T g x >P] 



(G 2 ) 



where (G 2 ) is the gluon condensate, obtained by the prescription 



rid /~ia 



{Vita ~ VliP Vua) (G 2 ) . 



(47) 



(48) 



In order to make predictions, we have to relate Gh in (|47| ) and the various loop integrals to physical 
quantities like m x , and an — fDy^nB- 

It should be noted that there are apriori other terms than the one in (H). There is one possible 
term where the field M q occurring in Fig. 7 may instead be attached to the light quark line in 
diagram in Fig. 8 (right). However, this term will not give contributions to D° — ► K°K°. Moreover, 
there is apriori a term where the field A a attached in Fig. 7 is instead attached to the light quark 
line in Fig. 8 (right). This term is identically zero. 

In the language of chiral perturbation theory, the term (^7j) can be interpreted as a counterterm. 
To be more specific, the (divergent part of the) counterterm has the Lorentz and flavor structure of 
the second line of (|47|) and is multiplied with a (divergent) coefficient adjusted to cancel the loop 
divergences obtained in Sect. 3. 
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6 Results 



In our numerical calculation we use the values of g and / obtained within the same framework 
in p8| , |27], ^2], [32], [33|. The coupling g is extracted from existing experimental data on D* — ► Dtt 
and D* —> decays. This analysis includes chiral corrections at one loop order and yields 

hO. 04+0.05 



5 = 0.27; 



0.02-0.02 ' 



leaving the sign undetermined. The one loop chiral corrections reduce the pion 



decay constant from f w = 0.132 GeV to / = 0.120 GeV pq| . In order to obtain the an coupling, we 
use present experimental data on D s leptonic decays. Namely, at the tree level there is a relation 
Id = /d s = olhI 'y/rnE- This relation receives 10 — 20% chiral corrections [22], 
experimental branching ratio D s — > fxUn and the D s decay width [15] one gets fjj 



2||. From the 
0.23 ± 0.05 

GeV and taking into account chiral loop contributions, we find an = 0.23 ± 0.04 GeV 3//2 (see 
Appendix [B]). Note that in [|^] the A = has been used, while we use the strict MS prescription 
A = 1 as in [28]. We put everywhere fi = 1 GeV ~ A x . 

For the Wilson coefficients c a ,b of @ we use c A = 1.10 ± 0.05 and c B = -0.06 ± 0.12 
calculated at the scale [i = 1 GeV with the number of colors N c = 3. Within the framework 
of "new" or "generalized" factorization, where nonfactorizable effects are taken into account in a 



phenomenological way, one uses the "effective values" c 



1.26 and c 



eff 
B 



-0.47. However, in 



this paper we calculate nonfactorizable effects in terms of chiral loops and gluon condensates, and 
therefore we use the values of |2q| . Due to the suppression of cb in comparison with ca, we do not 
include terms proportional to eg. [] We present our numerical results for the nonzero amplitudes 
in Table @. 





Mi[xW-< GeV] 




-0.42 


M 2 


-0.31 


M 3 


-0.62 


M A 


0.28 - 2.Ui 


M 5 


-0.81 


M & 


-0.61 


M 7 


-0.99 


Ms 


0.92 




-2.56 - 2.44i 



Table 1: Table of the one chiral loop amplitudes (see Fig. ^J), where M — J2 n M is defined in (pp|). In the last line 
the sum of all amplitudes is presented. It can be compared with the experimental result |.Mexp| = 3.80 x 10 -7 GeV. 



The imaginary part of the amplitude comes from the -F4 graph, when the 7r's or the K's in 
the loops are on-shell. All other graphs contribute only to the real part of the amplitude. The 
imaginary part of the amplitude is scale and scheme independent within chiral perturbation theory. 
This amplitude is also obtained from unitarity, and is valid beyond the chiral loop expansion. We 
also mention that the rescattering contribution, considered in [[ 
the one we calculate from graphs on Fig. ||. 

In order to cancel divergences one has to construct counterterms. In our case, this is described 
at the end of section 5. Generally, one can do that by using the symmetry arguments, as it has 

Even if the "new factorization" values had been used, the cb part of weak interaction would be suppressed by 
1/3 compared to the ca one. 



1C] is the same contribution as 
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been done in [27, 28] for the semileptonic decays of heavy mesons and D* decays. In the case of 
D* [g8| it was estimated that the contribution of counterterms is not substantial. 

To obtain the D° —* K°K° amplitude due to gluon condensate we have to know the coupling 
Gh- In addition, we have to find the tensors S in (42, ^3|) (and thereby the T"s) for K°K° in the 



final state. We find to lowest order for the parts of T g ^ p : 

1 2 

s % = - S » = -J2( ms - m d) (P + P)ti ; Sff = Y^{m s + m d ){p-p)^ , (49) 



where p and p are the momenta of K° and K° respectively. From (46) we see, that the momenta 
will be contracted with = p^/Md; where pn = p + p. It is important that and have a 
different momentum structure than (K°K°\ Vp\0) in Eq. (||), and they will give a nonfactorizable 
contribution to D° — ► K°K° proportional to (G 2 ), while Sjf does not. Note that T^ ,p of Eqs. 
(41,42,47) do not contribute. We find the gluon condensate contribution: 

M(D° -> K°K°) m = CA { Gml) {ms ~ ™ d) ^ 5 g f D (50) 



where: 



5 G = N C ^^ , G H ^f3^, B g = m^{I G1 -I G2 ) = ^-. (51) 

When we take into account the various relations between the loop integrals (Ps) and Gh, we find 
that /3 ~ 1/4. Using the values @ (|G 2 ) ~ (333 MeV) 4 , m x = 200 MeV, and m s ~ 150 MeV, 
we obtain the numerical value: 

M{D° -» K°^ ) (G 2) ~ 0.87 x 10" 7 GeV ; (52) 

which is also of the same order of magnitude as the experimental value. 

Adding both the chiral loops and the gluon condensate contributions, we obtain the total 
amplitude to C(p 3 ) 

M Th = (-l.7-2.4i) x 10" 7 GeV , (53) 

or in terms of branching ratio 

B(D° -► K°K°) Th = (4.3 ± 1.4) x 10" 4 , (54) 

where the estimated uncertainties reflect the uncertainties in the input parameters, especially an- 
Around the charm mesons mass region there are many resonances. One might think that their 
contribution will appear in this decay mode, either as scalar resonance exchange like in [|| or as 
K* exchanges |, @, Within our framework they would appear as the next order contribution 
(C(p 5 )) in the chiral expansion . This is, however, beyond the present scope of our investigations. 
It is interesting to point out that the effects we calculate, both from chiral loops and from the gluon 
condensate, are results of the SU{3) flavor symmetry breaking. In the limit of exact symmetry 
both contributions will disappear. 

We can summarize that we indicate the leading 0(p 3 ) nonfactorizable contributions to D° — > 
K°K°. The chiral loops and gluon condensates give the contributions of the same order of magni- 
tude as the amplitude extracted from the experimental result. 
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A List of integrals from chiral loops 



Here we list the dimensionally regularized integrals needed for evaluation of xPT and HQEFT 
one-loop graphs shown on Fig. f|: 



d 



4-e, 



1 



1 



(27r) 4 - £ q 2 - m 2 



ifi 



d*~ € q 



1 



16vr 2 
1 1 



h(m), 



(2vr) 4 - e (q 2 - m 2 )(q ■ v - A) 16vr 2 A 



with 



(771^ \ — 
— —J — m 2 A 



m 



J 2 (m,A) = -2A 2 ln(- T ) - 4A F[ — ) + 2A (1 + A), 



(55) 
(56) 

(57) 
(58) 



where A 



7 + ln(47r) + 1 (in calculation A = 1), while F(x) is the function calculated by 



Stewart in [28], valid for negative and positive values of the argument 



X 



X 

Vx^T 



TT 



tan 



In (^x + \/ x 2 — 



\x\ < 1 
b| > 1 . 



(59) 



The other integrals needed are 



(27r) 4_<E (q 2 - m 2 )(q ■ v - A) IQtt 2 



[I 2 (m, A) + h(m)], 



(60) 



with 



<*7 (2^ (, 2 ~ J)t v-A) =T^ A ^ ^ + ^ *W ' < 61 > 



JxCm, A) = (-m 2 + ^A 2 ) In + 1(A 2 - m 2 )F (J) 



^A 2 (l + A) + l -m 2 (2 + 3A) + 2 -m 2 - ^A 2 



J 2 (m, A) = (2m 2 - % -A 2 ) In - ^(4A 2 - m 2 )F 



AJ 



+ ^A 2 (l + A)-^m 2 (l + 3A)-^m 2 + lA 2 



(62a) 



(62b) 



The functions Ji(m, A), J^im, A) differ from the ones in Boyd - Grinstein list of integrals p^] by 
the last two terms in ( |62| ) that are of the order of 0(m 2 , A 2 ). These additional finite terms originate 
from the fact that 77^" is 4 — e dimensional metric tensor. 
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The chiral loop integrals needed are 
d 4 - e q 1 



m 

i/j, e 



1 



(2vr) 4 - e ((q + k) 2 - m 2 )(q 2 - m 2 ) 



d 



4-e, 



16vr 2 



(2vr) 4 - e ((q + k) 2 - m 2 ){q 2 - m 2 ) 16vr 2 
d A ~ e q q^q u _ Wk v 

(2^) 4 ~ e ({q + k) 2 - m 2 )(q 2 - m 2 ) ~ "16V 2 



N (m,k 2 ), 
iVi(m,A; 2 ) 



1 jfe*» 



2 16vr 2 



AT (m,A; 2 ), 



N 2 (m,k 2 



167T- 



■N 3 (m,k 2 ), 



where 



iVo(m,A; 2 ) 
iY 2 (?n,A; 2 ) 



-A + l-tf — +ln 



m 



ctG^ ^ J sign(/i ) , 



■■ 2 ^ - 1 ^ 2 — ) A - i{ I(8m 2 + k 2 ) [l - iff f4) 
6/ 2l3 v ; 2 W/ 



7V 3 (m,fc 2 ) = - m 



—m 
3 
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k 2 + (m 2 + 



A- 



In 



A 4 



i7r0 (~^2") sig 11 ^ 2 ))} > 



and 



2(1- y^/a - larct 



an 



2 1- -y/l-A/l 



In 



VI - 4/a + 1 



yr^i/a-i 



mQ(a - 4) 



< a < 4 



otherwise 



while m 2 is assumed to be positive. 



(63) 
(64) 
(65) 

(66) 
(67) 

(68) 



(69) 



B D meson decay constant 

Here we list results for onedoop chiral corrections to D meson decay constants and use them to 
obtain coupling an from experimental data. The onedoop chiral corrections have been calculated 



in [p2| , |27H using A = 0, while the leading logs have been obtained already in [32], |34|] 

;L i + 3" 

327r 2 / 2 (j2 Il(ym ^ + ^(ra^) + \ Il ^ mr i) 



fD 



1 + 32n 2 f 2 {^ C ( A D*D,m n ) + C(A D * D , run) + ^C(A D * D ,m V/ 



an 



~32^( 2/l(mx) + ^ /l(mr ' ; 



(70a) 



(70b) 



16 



where C(A, m) = Ji(m, A) + A^Ji(m, A), while Ji(m, A) and Ii(m) can be found in appendix 
0. Using / = 120MeV, /x = lGeV and A = 1 one gets the numerical values 

aH ^1 + 0.18 - 0.37c; 2 ) , (71a) 




,1 + 0.35 + 0.38# 2 ) . (71b) 

To obtain the an coupling we use experimental data on decays of D mesons into leptons. From 
the experimental value for branching ratio B(D S — > [iv^) = (4.6 ± 1.9)10 -3 and the D s decay time 
TDs = (0.496 IJSJS) • 10" 12 s one gets f Ds = 0.23 ± 0.05 GeV. Using this value and g = 0.27 |f in 
(IfTbD we get a H = 0.23 ± 0.04 GeV 3/2 . 

Using a H = 0.23 ± 0.04 GeV 3 / 2 in (|71a| ) and # = 0.27 we can also calculate f D = 0.194 ± 0.045 
GeV, where the uncertainties are due to the uncertainties in an- The average value for the ratio 



Ids/Id = 1-19 is in fair agreement with the recent lattice results [29|. 



C Heavy-light quark loop integrals 



The integrals entering heavy quark loops like the ones in Fig. 8 are of the form: 

d d k 1 1 



(72) 



(2ir) d (v.k)P (k 2 - m 2 )i ' 
Performing a shift of momentum integration combined with Feynman parameterization, we obtain 



R 



p.q 



l) 



(73) 



where 



K{n, r) = I <1\ 
o 



d d l 



y 



(2Tr) d (i 2 -m 2 - x 2 y 



One should notice that to obtain the result in (13), we have to do the identification 



8iN c G 2 H I HH = 1 , 



(74) 



(75) 



where Ihh is a logarithmically divergent loop integral given below. (There is also a similar relation 
for g.) One should notice that some authors use an extra factor ran, the mass of the heavy meson, 
in front of the right hand side of (13). Choosing the normalization in (|l3|), it means that a factor 
^rriH is included in the heavy meson field H v . For the left handed current in (18) and (|H]) we find 
that we have to identify : 



a H = -4iN c G D I H w , 
where Ihw is a quadratically divergent loop integral. 



(76) 



The regularization can be done in various ways (various cut-off prescriptions or by MS) and each 
regularization correspond to slightly different versions of this type model [ 1^, 18 , [L9|, |20|, ^T], Q . For 
instance, in the version we use, when soft gluon emission is included in (^5) above, gluon condensate 
contributions should also be included in loop integrals Ihh and Ihw, as it is for /„- in the light 
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sector [12, However, we will not go into these details here. Still, as in the pure light quark 
case, one obtains numbers in the right ball park by parameterizing the quadratic divergent integral 
as A 2 , and the logarithmically divergent integral as log(A 2 c /m 2 ). Anyway, using the expression for 
f n obtained in the xQM, we obtain from (|75|) to leading order 



G H 



2, /to 



x 



(77) 



It can be seen from Ward identities for the loop diagrams for, say Fig. 8 (left) that the quadratic 
divergence in Ihw is related to the quark condensate of the light quark, which is also quadratically 
divergent. Then, similar to (f77l), we obtain from ([7(1) to leading order 



G H 



Combining (f77|) and ( |78| ) we obtain 



a H 



2 m x a H 

m 
m 



(78) 



(79) 



which for the values m x = 200 MeV, f n = 131 MeV and (qq) = (-240MeV) 3 gives the value for an 
cited in Appendix B. Furthermore, using ( f77|) and (|78|) we obtain 



-2^4- ~ 0.25, 

(qq) 



(80) 



to be used in ( pOP and ( |5l| ) . 

Within dimensional regularization, the expressions for some vales of n and r are listed below: 



K(2,l) 
K(3,l) 
K(3,2) 



r(i - d/2) 



2(4ir) d/2 (to 2 ) 1 "^ 2 ' 
i r(2-d/2) 



4(4vr 



d/2 



> 2 )2-d/2 



T(3/2 -d/2) 



16(An) d/2 - 1/2 (m 2 ) 3 / 2 - d / 2 



where to = m x . From the properties of the T function it is easy to see that: 



K(2,0) = -AK(3,2) 



m 2 K(3,0) = (3-d)K(3,2) , 



(81) 
(82) 
(83) 

(84) 



Comparing with a cut-off regularization, we see that K(2, 1) is quadratically and K(3, 1) is loga- 
rithmically divergent. In a primitive cut-off regularization K(2,0) and K(3,2) appear as linearly 
divergent while they here appear as finite! 

Note also that some of the integrals (74) can be obtained as the limits of integrals listed in 
appendix ^ if one lets A — ► 0. Thus one has the relations 



if (2,0) 
if (2,1) 
if (2, 2) 



1 



32vr 2 A^o A 



lim —l2(m, A) 



k lim 

32vr 2 A^o 

i 



I 2 (to,A) + /i(to) 



32vr 2 



h{m) 



32tt 2 A^o 



lim A J2(to, A) . 



(85) 
(86) 
(87) 
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The loop integral's (the J's) are defined as: 



I HH = mK{3, 1) + K(3, 2) , 

I HW = K(2,l)+mK(2,0) , 

I G1 = K(3,l)+mK(3,0) , 

IG2 = K(3, 1) . 
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